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A THEORY FOR THE LATERAL RESPONSE OF AIRPLANES 
TO RANDOM ATMOSPHERIC TURBULENCE 
By John M. Eggleston 


SUMMARY 


The lateral motions of an airplane flying through continuous random 
isotropic turbulence have been derived in terms of (l) the transfer 
functions relating the motion in the various degrees of freedom to the 
yawing moment, rolling moment, and side force, (2) the statistical forces 
and moments at the center of gravity due to gust velocities acting on the 
lifting surfaces of the airplane, and ( 3 ) the power spectra of the three 
orthogonal components of gust velocity acting on the airplane along the 
flight path. The method takes into account the random variations of 
gust velocity across the span and along the fuselage. Solutions are 
given in the form of equations relating the power spectra of the angular 
motions of the airplane to the power spectra of the gust velocities. 

Three airplanes of different size are used to demonstrate the method, 
illustrate characteristic trends, and exhibit some simplifications pos- 
sible in the calculations. For these airplanes the effects of horizontal 
gusts (that is, gusts parallel to the flight path) and side forces due 
to gusts on the airplanes were found to be negligible. 

By using one of the example airplanes, a comparison is drawn between 
the present theory and several less comprehensive theories for calcu- 
lating the effect of gusts on wings of finite span and the effect of 
gusts on the motions of the complete airplane. 


INTRODUCTION 


The classical theory of stability and control of airplanes has been 
applied to the calculation of response to controls and response to gusts. 
In the calculation of response to gusts, however, the angle-of-attack 
distributions which produced the forces and moments applied to the air- 
plane by the gust have customarily been assumed to be equivalent to those 
resulting from a rigid-body motion of the airplane. This method was 
applied in NACA Report 1 and other early reports (ref. 1) in the study 
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of the longitudinal response of airplanes to gusts and in reference 2 
to the calculation of the lateral response to gusts. Because of the 
method of accounting for the effect of the gusts, these analyses were 
restricted to the consideration of the effects of isolated gusts of some 
specified shape. Most early developments in the theory of response to 
gusts were attempts to account for factors which were important in the 
calculation of loads. For example, the calculation of loads due to 
vertical gusts was extended in reference 3 to include the effects of 
unsteady lift and of flexibility. 


In the early methods of calculating gust response, the analysis did 
not include such factors as the penetration effect (a result of the 
penetration of a given gust at different times by different parts of the 
airplane) and the variation of gust velocities along the fuselage and 
spanwise directions. In recent years, some attempts have been made to 
include these effects without changing the basic method of approach. In 
reference 4 an approximate method was introduced for estimating the 
penetration effect in the calculation of longitudinal response to gusts 
by taking into account first-order effects of the time lag between the 
penetration of the gust by the wing and by the tail. The method used 
was similar to that introduced in reference 5 for calculating the effect 
of lag of downwash. A similar application of this concept to the calcu- 
lation of lateral response to gusts was presented in reference 6. 

A more realistic treatment of random disturbances in the atmosphere 
was made possible with the introduction of a method of calculating the 
statistical response of an airplane in terms of the statistical proper- 
ties of atmospheric gust velocities. This method, known as statistical 
dynamics, was first applied in reference 7 to the measurement of the 
longitudinal response of an airplane to random vertical gusts and has 
since been used extensively in longitudinal gust loads studies. 


A theoretical analysis for the statistical lateral response of an 
airplane to gusts was given in reference 8 where, in order to approximate 
the effects of the span, the turbulence was considered to be represented 
by a combination of rolling gusts and side gusts having arbitrary distri- 
butions along the flight path. With this representation of turbulence 
however, the question arose as to the relative magnitudes of the rolling 
and side gusts which would satisfactorily represent atmospheric turbu- 
lence. In reference 8 the assumption was made that an arbitrary gust 
distribution across the span could be represented by a constant anti- 
symmetric gradient which was assumed to be equal to the gradients existing 
m the corresponding distribution of side gusts. Experimental flight 
measurements of reference 9 indicated, however, that the measured gradi- 
ents of vertical gusts across the span did not agree with the theoretical 
calculations of reference 8, and the assumption of nonisotropic atmos- 
p eric turbulence was employed in order to fit the theory to the flight- 
test results. ^ 
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A more realistic method of utilizing statistical dynamics to take 
into account random gusts both along the flight path and across the span 
was reported in reference 10 in outlining a means of calculating the 
longitudinal response of an airplane. The analysis was consistent with 
isotropic turbulence and required only the statistical gust distributions 
measured at one point in the turbulence. 

Therefore, with the advent of statistical dynamics and more realistic 
means of analyzing the random distributions over the airplane as well as 
along the flight path, the original methods of treating gusts no longer 
appear to be justified. Furthermore, although it would be desirable to 
continue to treat the effects of gusts as equivalent to certain rigid- 
body motions of the airplane, such treatment is not consistent with 
statistical gust inputs. The forces and moments on the airplane rather 
than the resulting motions are statistically related to the gust input; 
that is, a random distribution of gusts across the span and along the 
fuselage of an airplane will produce changes in forces and moments which 
are related to the gust velocities along the flight path only in a sta- 
tistical manner. The purpose of this paper is to establish the statis- 
tical relations between gust velocities and the forces and moments 
affecting the lateral motions of the airplane and to use these forces 
and moments to compute the lateral motion. 

In preparation for the calculations of the present report, two pre- 
ceding papers have been published. One of these papers (ref. 11) con- 
tains a method for calculating the side force and yawing moment on the 
fuselage and tail of an airplane in continuous gusts. The other (ref. 12) 
presents the lateral forces and moments on a wing due to three-dimensional 
random turbulence for a number of wing load distributions and for wings 
of various spans. This work (ref. 12) is an extension of the theory 
given in reference 10 for defining the lift on a wing due to random 
turbulence . 

This paper presents a method of utilizing the statistical forces and 
moments derived in references 11 and 12 for the calculation of the lateral 
motion of an airplane in any lateral degree of freedom due to atmospheric 
turbulence . First, the governing equations are derived and the power 
spectral solutions of the motions are given in terms of linear airplane 
transfer functions, statistical force and moment coefficients, and power 
spectra of the three orthogonal components of atmospheric gust velocity. 
General expressions are given or referenced for each of these independent 
functions. Next, solutions for three airplanes are shown with trends and 
possible simplifications noted. Finally, a comparison is made between 
the results of the present theory and those of earlier theories, both by 
comparing the gust distributions required to produce equivalent airplane 
motion and by comparing the airplane motions for the same assumed gust 
characteristics . 
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SYMBOIS 

A constant used in references 6, 8, and 9 to indicate relative 

amplitudes of power spectra; aspect ratio (table II) 

B matrix of equations of motion of airplane 

b wing span 

D nondimens ional operator, ^ ^ 

E(wg) plot of °f reference 12 

f arbitrary function of time 

F Fourier transform of a quantity 

hp pressure altitude 

h z height of center of pressure of vertical tail above X 

stability axis of airplane 

i-\R 
k n = cox n /u 
k' = cdL/U 

arbitrary constant 

nondimens ional radius of gyration about X stability axis 
nondimens ional radius of gyration about Z stability axis 
nondimens ional product of inertia 

tail length measured from center of gravity to center of 
pressure of vertical tail 

integral scale of turbulence 

arbitrary terms which may take on values 1, 2, 3 • • • 
rolling velocity, (f) 


K 

K x 

K z 

K X Z 

h 

L 

m,n 

P 


1 


dynamic pressure 


NACA TN 3954 


5 


r yawing velocity, f 

R real part of complex quantity 

S wing area 

s profile height (see eq. (44)) 

t time 

T particular value of time 

U relative velocity between airplane and general air mass 

u velocity along X stability axis 

v velocity along Y stability axis 

w velocity along Z stability axis 

W weight of airplane 

X,Y,Z three orthogonal stability axes (see fig. l) 
x,y,z coordinate position with reference to stability axes 
Cj^ lift coefficient, 

qS 

Cj rolling-moment coefficient, momen t 

qSb 

C n yawing-moment coefficient, Yavin g moment 

qSb 

Cy side-force coefficient, Sfde force 

qS 

C Do coefficient of drag at zero lift 

a angle of attack 

P angle of sideslip 

3' = b/L 

r effective dihedral angle 
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0 

A 

P 


P 

a 

* 


CJD 


0 )' 

“'n 

%lf 2 


flight-path angle 
angle of pitch 
taper ratio 

relative density factor, j^ s 

density of atmosphere 

sidewash angle 

angle of roll 

angle of yaw 

circular frequency 

reduced frequency, aib/U 

natural reduced frequency 

power spectral density of a function f 

cross-power spectral density of functions f]_ and f 2 


C.P. cross power 

Stability derivatives of airplane are indicated by subscript notation; 
for example , 



Subscripts : 

W wing 

F fuselage 

T vertical tail 

FT fuselage and vertical tail 


g 


gust 
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0 trim condition 

j,k,m generalized matrix subscripts 

Superscripts : 

* conjugate of a complex number 

Bar over quantity denotes mean value. 

Dot over quantity denotes derivative with respect to time. 


Matrix notation: 


- 

i 


absolute value of a quantity (determinant of a matrix) 
rectangular matrix 
column matrix 


LATERAL EQUATIONS OF MOTION 


Systems of Axes and Description of Turbulence 


The motions of an airplane flying through continuous atmospheric 
turbulence are defined as angular displacements, rates, or accelerations 
about a set of stability axes (with respect to the general air mass, the 
X stability axis is always alined with the projection of the resultant 
mean velocity on the plane of symmetry) whose origin is fixed at the 
center of gravity. These stability axes have some mean orientation with 
respect to the body sixes of the airplane denoted by the angle a Q in 
the XZ-plane. The general orientation of the stability axes of the air- 
plane with respect to the gust velocities and with respect to the body 
and earth axes is shown in figure 1. 


At each position along the flight path of the airplane the instan- 
taneous velocities of the local air mass, hereinafter designated as gusts, 
are instantaneously defined by three orthogonal components of velocity 
alined with the stability sixes of the airplane at that position. The 
gusts are assumed to be essentially isotropic in nature, and their statis- 
tical characteristics remain unchanged regardless of the motion or 
orientation of the airplane. This assumption is implicit throughout 
this paper. 

Furthermore, the assumption is made that, during the small interval of 
time in which each gust acts on the airplane, the gust velocity does not 
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change appreciably. On this basis the three components of gust velocity - 
alined with, but of opposite sign to, displacements along the stability 
axes of the airplane - are functions only of their space position along 
the flight path. With reference to the center of gravity of the airplane, 
these gusts are denoted by 


u g = ug(x,y,z) 
v g = v g (x,y,z) 

v g = w g( x >y' z ) 

The orientation of the gust axes is shown in figure l(a) . 

The assumption of isotropy further relates these velocities only to 
the extent of 



where the bar denotes the mean value. 


System of Equations 

The lateral equations of motion of an airplane defined in the 
stability system of axes are generally well known. Based on these equa- 
tions of motion, the frequency response of an airplane in its three 
lateral degrees of freedom is defined by the matrix 


\ 

0(a)) 



Cn 

Cy 


\ 

Cj(od) 


> = 


fw> 

£(<*>) 

c Y 

s 

C n (oo) ? 

P(a>) 




fM 

Gy 


Cy(co) 

^ 4 


In equation (l) the ratios 0/Cj, \|r/C n , and so forth, are complex transfer 

functions of the airplane and relate the response of the airplane in 0, 
i| |r, and p to a sinusoidal rolling moment, yawing moment, or side force 
of unit amplitude. These moments and side force are functions of the gust 
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velocities referenced to a point (the center of gravity of the airplane) . 
The necessary phase relations between these forces and moments and the 
gust velocities at the center of gravity are taken into account by 
expressing C^, C n , and Cy as complex quantities. In terms of the 
gust velocities at the center of gravity, the lateral-force and moment 
coefficients may be expressed as follows: 


C j(cd) 


<C n (o>)> = 


^(a>) 

w j 


u 


g 


'g 


g 


Sl(») 

U g V g g 


u g (o>) 


<Vg(co)^ 


Cy(co) 

^ J 


u s 


V) 

V g 


^(o>) 

Wg 


Wg(o)) 

V 


( 2 ) 


In equation (2) the elements of the rectangular matrix are transfer 
functions relating the forces and moments to the sinusoidal components 
of gust velocities of unit amplitude measured at the center of gravity. 
Inasmuch as some of these forces and moments arise as a result of gust 
velocities distributed over the airplane, which are, in turn, related to 
the gust components at the center of gravity only in a statistical sense, 
these transfer functions are likewise statistical in nature and are later 
expressed in the form of power spectra. 


With the foregoing restrictions the complete relationship between 
the frequency response of the airplane in the three lateral motions and 
the gust velocities is given by the matrix equation 


0 (®) 

3 (ao) 


£ 

cz 

± 

Cz 


Cz 


1 

Cn 

JL 

Cn 




£ 

Cy 

JL 

cy 

_P 

Cy 


zl 

u g 


■'ll 

g 


U, 


U 


g 



“I 


r 

£ 

£ 


u g (o>) 

v g 

w g 


£n 

Cn 

< 

-A. 

bO 

> 

v e 

w g 


Cy 

V g 

£y 

w g 


Wg(o)) 


(3) 


Thus, the moments due to the gusts and those due to the resulting 
motions are superposed. The equations of motion of the airplane are 
assumed to he linear within the range of this resulting motion. 
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Power Spectral Response 

Inasmuch as turbulence in the atmosphere appears to he random, the 
statistical response of an airplane may, in general, he predicted in 
terms of the power spectral densities and mean-square values of the gusts 
in the atmosphere. An excellent summary of this field of dynamics is 
presented in reference 13 . From the relationship (eq. ( 3 )) between the 
gust velocities and the motions of the airplane in the frequency domain, 
the relationship between these quantities in terms of power spectra and 
cross-power spectra may be readily defined. 

A quantity (herein denoted, in general, by the symbol f(t)) which 
has a random variation with time in the range -T ^ t S T and is zero 
elsewhere is related to its Fourier transform F(o)) by 



where, in turn, 

f(t) = | r F(a>) e ia)t da) (5) 


The power spectral density of the quantity is defined (see, for 
example, ch. 6 of ref. 14 or appendix of ref. 9 ) "by 


if (a)) 


= lim — F*(cd) 

T-»oo T 


F(o>) 


= lim i | F(o>) | 2 
T->°° 1 


( 6 ) 


where asterisks are used to denote the complex conjugate of a quantity 

and the term — indicates that the spectra are symmetrical and only 
T 

positive values of a) are considered when taking the average value. 

When two related quantities are measured concurrently (for example, 
when the yawing and rolling moments due to each gust component are 
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measured at the same values of time) , a cross-power spectral density 
exists between these two quantities. Where the frequency content of 
these two quantities is defined by 



> 


(7) 


the cross-power spectral density is defined by 

®f _ f _ = !i m ^ F 1 *(oj)F 2 (a)) 
1 d T->oo 1 


or 


®f 9 f, = lim A F 2 *(a3)F 1 ((o) 

^ X T ->00 

= 'If ^2* (8) 

The significance of this cross power lies in the fact that, if the quan- 
tities are related, a certain interchange of energy takes place between 
them and the amount of this interchange is defined by the cross -power 
spectrum. 

In the special case where F]_ and F 2 are the same function and 

are not removed in time or space (that is, f^(t) = f 2 (t)), equations (6) 

and (8) are identical and the cross power of f with itself is, in 
fact, the power spectrum. 

These definitions may be found in most texts and in a number of 
papers on power spectral analysis, among them references 9 , 13 , and 14. 

On the basis of the preceding discussion, it may be shown that if each 
complex equation of the system defined by equation ( 3 ) is multiplied by 
its conjugate equation, divided through by the quantity T, and the 
limit as T — taken on both sides of each equation, then the system 
of equations, by definition, is in the power spectral form. The power 
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spectra and cross-power spectra of the variables may be recognized and so 
denoted, the spectra and cross spectra of the gust velocities being 
unknown variables. An example of this operation is shown in appendix A. 


By using this or any other suitable method of power spectral analy- 
sis, the system of equation (l) may be defined in terms of power and 
cross-power spectral densities such that 


<$• 




= 


0 

2 

9* 

2 

1 

c z 


C n 


C Y 

♦ 

2 

* 

2 

♦ 

Cz 


Cn 


CY 

3 

2 

J_ 

2 

3 

c z 


C n 


C Y 


$Ci 




<$ P + 2R 

u n 


\*fi (j\*t 


\?l C n 


\ C n/ C Y 


\c Y c z 


(±\*±_ 


C i C n \C n Of \Cy/ C i 


3 \* 3 / 3 \* 3 / 3 \* 3 


?l C n \ C n C Y v Y/ C Z 


®C Z C n 




Vi 


(9) 


where 2R means twice the real part of the product of the two matrices. 


On the right-hand side of equation (9) , the product of the first 
two matrices is the contribution of the power spectra of the forces and 
moments, and the product of the second two matrices is the contribution 
of the cross-power spectra between the moments and forces produced by 
any one gust component and also the cross powers between the forces and 
moments produced by any two gust components. The powers and cross powers 
of the forces and moments may, in turn, be related to the powers and cross 
powers of the gust velocities. In so doing, a great simplification is 
obtained when atmospheric turbulence is considered to be isotropic 
within the frequency range affecting the response of the airplane. 


By virtue of the assumption of isotropic homogeneous turbulence, 
the phasing between the u g and w g components and the v g and 
w g components is unrelated and the cross power of the components is 
zero. The phasing between the two gust components lying in the XY -plane 
of the airplane, that is, u g and v g , is unrelated when these veloci- 
ties are measured at a given spanwise station (for example, along the 
airplane center line) but is related when u g at one spanwise station 
is compared with v g at another spanwise station. However, the effects 
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of the cross power between these two components are assumed to be 
insignificant, and justification for this assumption is made in a sub- 
sequent section. Thus, 


$> u v (to) = 
18 


#v w (to) = 
88 


K u (“) = 
88 


( 10 ) 


By the method of multiplying by the conjugate and taking the limits, 
the matrix relationship of equation (2) is defined in the form of power 
spectra by 






<h 

2 

C l 

2 

C l 

2 

u g 


V g 


W g 


C n 

2 

c n 

2 

Si 

2 

u g 


v g 


w g 



( 11 ) 



Cy 

2 

c Y 

2 

c Y 

U g 


v g 


w g 



Likewise by cross -multiplying the equations in the matrix of equa- 
tion ( 2 ) to conform with the definition of a cross-power spectrum, the 
cross spectra between the forces and moments are given by the matrix 



A proof of the relationships of equations (ll) and (12) is given in 
appendix B. 
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Correlation Between Airplane Components 

Generally, a force or moment acting on or about the center of 
gravity of the airplane is due to the penetration of the gust by each 
component of the airplane, and the correlation or phasing between the 
development of these forces and moments must be included. Complete 
generalization to all possible combinations of components would lead 
to lengthy and unnecessary expansions; however, all lifting surfaces 
which might reasonably be expected to contribute significantly to the 
lateral motion of the airplane as well as the components of the gust 
velocity that will significantly influence the forces and moments on 
these lifting surfaces are incorporated. For some unusual configura- 
tion where other lifting surfaces are thought to be important, the 
method of including these factors is clearly indicated. 

The force and moment coefficients which include only those func- 
tions deemed pertinent to the lateral motion of the airplane are defined 
by 

Cz(ug,vg,w g ) = (Ci)w(ug,vg,w g ) + (Cz)T(vg) 

C n (ug,vg,w g ) = (C n ) W (ug,w g ) + (Cn)FT(vg) ’ 

CY(ug,Vg,w g ) = (C Y )pr ( v g) 


In terms of the complex force and moment derivatives, the force and 
moment matrix is thus defined by 



( 14 ) 


Y 
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By the same method as was applied to equation (2) to obtain equa- 
tions (ll) and (12) (multiplying by the conjugate and taking the 
limits), the power spectra and cross-power spectra of the force and 
moment relationships of equation ( 14) are immediately definable as 


/ N 

_ 


1 

1 

*( c i)t 


^( C n)w 

— ».r~ 

II 

^ (Cn) FT 


$ ( c y) ft 
\ > 



u. 


g 


W 


0 


'n 


u 


■g 


W 


£i 

v g 

Cz 


'g 


Cn 


FT 

2 

FT 


w, 


g 


0 


w. 


0 


*U 


< 


g 


g 


w. 


g 


( 15 ) 


®(Cz)w(Cz)t 


0 

©W©T 

0 




(Si\*(9r) 

WM u g/w 

0 

WN 

WMvw 


®(Cz)w(C n )FT 


0 

(Sifp) 

\ V g/W\ v g/FT 

0 


®(cz)w(Cy)ft 


0 

WM v g/JT 

0 


®(Cz)T(Cn)w 


0 

0 

0 

% 

< ®(Cz)T(Cn)FT 

> = 

0 

ill 

0 


*(Cz)t(Cy)ft 


6 

ill 

0 

>j 

^(CjwCCnJFT 


0 

0 

0 


*( c n)w( c Y)FT 


0 

0 

0 


®(C n )FT(CY)FT 


0 

K (-) 

\ v g/FT\ v g/FT 

0 



(16) 
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The complete power-spectral-density relationships including the 
lateral moments and forces due to the various lifting surfaces of the 
airplane are similar to those given in equation (9) with the quantities 
defined in equations (15) and (l6) inserted. When the left-hand side 
of equation ( 15 ) is denoted hy {^Cj} and the left-hand side of equa- 
tion (l6) , by j<DCjC k } the complete relationship is defined by 


IW = 


0i 


Cj 


P Cj 


+ 2R 


Vi\*0i 

,Cj c k 


j^k 


where 0]_ is (j) } (f ) 2 is \|r, and so forth. The complex force and 

moment derivatives contained in equation (l4) were derived in such a 
manner that the phase relationships between the moments and the gust 
disturbances are given with reference to gusts acting at the center of 

fe 

c j 

in equation (l). As an illustration, the total response of the airplane in 
roll is obtained in the following form: 


gravity. For this reason, the transfer functions 


axe those given 


4 


0 

2 

0 

57 

^ ( C ^ ) W + ^^TJ + 

c Y 

$ (c y )ft + 

C n 


^(C n )w + ^(C^FTf 


2R 


'0\*0 /0\*0 

Cl) Cl ®( g z)w(Cz)t + C n $ ( c z)w(C n )w + 


\Pi ) c n °(C2)w(C n )FT + (c j) Cy <I> (Cz)w(Cy)fT + 


! / (f) / (h \ ^ 0 

~l) ^ $ (Cz)T(C n )FT + (c^) °( c z)t(Cy)fT + 


(I;) ^ *< c n)fT(CY)ra + 3 zer0 teras 


(17) 


Other degrees of freedom may be obtained directly by replacing 0 with 
the appropriate symbol, such as \|r or 3 . 
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By substituting in equation (17) the expressions for the power and 
cross-power terms of the force and moments from equations (15) and (l6), 
the motion is defined in terms of the gust velocities: 


JL 

2 

<t„ + 

£ 

2 

<k r + 

£ 

u g 

u g 

v g 

V g 

w g 


■v, 


g 


(18) 


where 


£ 

u g 


0 2 

Cl 


8c 1 


5u, 


■g 


W 


^ii 


8Cn 


8u, 


g 


+ 2R 


W 


8Cz\*/8Cn\ /£\*0_ 
8u g Jw\8 u gy W \ c 2/ c n 


(19) 


0 

2 

0 

2 

8c* 

2 

| 

0 

2 

8c 1 

2 

1 

0 

2 

dCn 

2 

_L 

0 

2 

8Cy 

v g 


Cl 


8v g 

w 

Cl 


8v g 

*r 

T 

Cn 


8v g 

1 

FT 

C Y 


8vg 


£ 2 + (££* (££\ (£\- 
]\8 v g/w\8vg/ T c z \8v g / FT \8v g / FT \c n / Cj 


2 

FT 


+ /£lYY£ n\ 

r g/wWg/FT \8vgyT\8vgJpr 

* 




r g /FT 7 v g /FT vn ) c Ti 

A*£ 


( 8c A*/ 8Cy \ + f 8Cz\*/8Cy\ 

V^ v g/W\8v g /FT \8vg/T\8v gy t T 


\Ci) c n 

C Y 


(20) 


0 

2 

0 

2 

8c z 

2 

J. 

0 

2 

Sc n 

2 

+ 2R 

/8cA*/8c n \ /0\*0 

w g 


C l 


8w g 

T 

w 

c n 


8w g 

W 

\8 w g/w\8 w g/w\ C Z j C n 


( 21 ) 


Thus, the response of the airplane in any degree of freedom depends 
on the power spectra of the three gust components, the eight complex 
coefficients of equation (14) which relate the gust velocities to the 
forces and moments, and, for each degree of freedom, the three complex 
transfer functions relating the motion to a sinusoidal rolling moment, 
yawing moment, and side force of unit amplitude. These coefficients, 
transfer functions, and power spectra are defined in subsequent sections 
of this paper. 


18 


NACA TN 3954 


Transfer Functions 

The lateral equations of motion usually accepted for the rigid 
dynamics of an airframe for which noncoupling with the longitudinal mode 
and all assumptions associated with small perturbations about a trim 
condition are valid are given by the matrix 


where 




r 

0 


Cj(ug,Vg,Wg) 

- * 

> = 4 

C n (ug,v g ,wg) > 

3 


CY( u g;Vg,w g ) 
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(22) 


2^ D - C Y)3 

(23) 


The solution, in terms of the rolling moment, yawing moment, and side 
force may be found by methods of operational calculus. One method is^ 
used by Mokrzycki in reference 15 where, except for the sign of K^, 

the identical equations of motion are used but with moments from 
aerodynamic-control surfaces as forcing functions. Another means of 
solution is through matrix algebra where the solution as given by equa- 
tion (l) is of the form 


0 

r-n-1 
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C l 
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The prime denotes the transpose of the matrix B. Such a method is 
found in reference 1 6 and other texts on matrix algebra. The solution 
by either method gives the complex transfer functions with real constant 
coefficients of the powers of D. With the common denominator denoted 
by | B | (equivalent to E& of ref. 15) , the transfer functions are 
given by 


£ 
c Z 


VVd 3 - + %)d 2 ♦ (| ♦ 2 ^ - | C^Cj^D - 


c np C L tan 7 
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IB 
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In addition to the angular-displacement transfer functions given 
here, the response in terms of the derivatives of the angular displace- 
ments may be obtained by multiplying both sides of equations (26) to (3*0 
by the operator D. For example. 


For the frequency response of these transfer functions, the operator 


where a) 1 is the reduced frequency in terms of radians per span length. 
Thus, the transfer functions relating the dynamics of the airframe to 
the lateral forces and moments are obtained from the inertial character- 
istics and stability derivatives of the airframe under steady flow con- 
ditions (that is, nonturbulent flow conditions). 


The moment and force coefficients which relate the rolling moment, 
yawing moment, and side force to a sinusoidal gust velocity of unit 
maximum amplitude acting in planes of (but of opposite sign to) the 
reference axes of the airplane are, in general, frequency-dependent, 
and only their statistical average properties are definable. Their 
derivation is based on generalized harmonic analysis and basic aero- 
dynamic theory. Where the derivation is lengthy, only the results of 
referenced papers are given or indicated. Furthermore, only those eight 
coefficients necessary to equation ( lb ) are defined. All others are 
considered as having small or negligible influence on the results. 

Rolling moments .- Four contributions to the rolling moment on the 
airplane are considered. They are the rolling moment of the wing due 
to all three components of gust velocity and the rolling moment due to 
the side gust acting on the vertical tail. 

In reference 12 are derived the power spectral densities of the 
rolling moments of a wing in isotropic turbulence due to horizontal, 
vertical, and side gusts (that is, the three components of gusts con- 
sidered herein) . The theory considers random turbulence across the span 
as well as along the flight path, and the results are given for several 
span-loading distributions and for a range of values of the ratio of 


i = u n i = u 
C Y b C Y b 


2h(% 2 % + %z%)d 5 + sKS - ill 


( 36 ) 
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Complex Moment and Force Coefficients 
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span length to the integral scale of atmospheric turbulence. The power 
spectra of rolling moment given in reference 12 are not repeated herein 
but are simply related to the notation of this paper: 


Cl 

v g 


2 

W 


®u g (“) = ®C 1 


(due to horizontal gusts, ref. 12) 


(38) 


Cl 


u, 


g 


< (to) = 0 q 7 (due to vertical gusts, ref. 12) 

w § ^ 


(39) 


For the side gust Vg, the random variation along the flight path 
was assumed to be uniform across the span, an assumption that is partic- 
ularly valid for airplanes having wing loading due to sideslip concen- 
trated over a small region of the span near - the plane of symmetry. Such 
an assumption is likewise made herein, so that 



is assumed to be independent of frequency. 


(40) 


The vertical-tail contribution to the rolling moment 


is defined by 


(Cj)tM 


\ dp/T qSb 
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1 + 8aV"i(W u ) 
dp I 


and, hence, 



Y Pj T bU 




-i(Z t o>/u) 

e 


(41) 


Thus, the rolling moment due to side gusts may he directly calculated 
from the known characteristics of the airframe in still air, and the 
rolling moment due to horizontal and vertical gusts may he obtained in 
power spectral form from the plots of reference 12. The fact that the 
real and imaginary parts may not he ascertained for the rolling moment 
due to Ug and Wg components is shown to he irrelevant; that is, the 
real and imaginary parts, as such, are not necessary to the calculation 
of the cross -power terms involved in calculating the response of the 
airplane . 
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Yawing moments .- The yawing moments due to the u g and w g compo- 
nents of the gust are primarily due to the unsymmetrical distribution of 
these gust velocities across the span of the wing. Since this unsymmet- 
rical distribution is the same characteristic of turbulence that produces 
the rolling moment on the wing, the yawing moment may be fairly accurately 
related to the rolling moment through aerodynamic relations. As given 
in reference 12, the expressions for the yawing moment of the wing due 
to small horizontal- and vertical-gust disturbances are 



nr\ /C_z 


c Ww\ u s/ w 


(cd) 


( 42 ) 





(43) 


where C nr , C; p , and so forth, are stability derivatives of the wing 

only at the trim (mean) angle of attack of the airplane and, as such, are 
not functions of frequency. Any yawing moments of the wing due to side 
gusts are neglected. 


The yawing moment of the fuselage and vertical fin of an airplane 
due to penetration into side gusts Vg has been derived in reference 11 
by using a simple profile shape to define the lift distribution over the 
fuselage and vertical tail due to sinusoidal side-gust components. As 
given in reference 11, 
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where x Q , x 1 , x 2 , s Q , and S]_ are geometric dimensions of the 

profile shape and frequency appears in the form 

*n = — 

^ u 


(n = 0, 1, 2) 
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For any given flight condition for which the exact value of Cn^ 

is known, it is recommended in reference 11 that small adjustments in 
the profile shape he made so that at zero frequency 



(oJ=0) 


he- 


U \ n 3 JFT 


Side force.- The side force on an airplane in gusts is considered 
to he primarily due to the effect of side gusts on the fuselage and 
vertical fin (the profile shape) of the airplane. By using a simple 
profile shape and slender-hody theory to define the lift distribution 
over the fuselage and vertical tail due to sinusoidal side gusts, the 
side force per unit side gust as derived in reference 11 is 


VgyFT 


2rt 

2sq 2 

SU 

CVJ 

/ 
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1 - (1 - iko)e 
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-ike 



( 45 ) 


where the notation is defined in the preceding section. Again it is 
recommended that small adjustments he made in the representative profile 
shape so that at zero frequency 


3d (£0=0) 

i V g/FT 



FT 


where (C v \ is known from flight or wind-tunnel tests for the desired 
y 1 3jFT 

flight conditions. 


POWER SPECTRA OF GUST VELOCITY 


As discussed in reference 17, the relationship (correlation) between 
gust velocities measured at two points in isotropic turbulence may be 
completely defined in terms of their lateral and longitudinal components 
(components measured perpendicular and parallel, respectively, to the 
vector distance between two points) . As treated herein, the airplane 
senses only the lateral components of the Vg and wg gusts and both 
the lateral and longitudinal components of the Ug gusts. (See preceding 
section and ref. 12.) 
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Both theoretical methods (ref. 17) and flight measurements (ref. 18) 
have shown that the following expression approximates the spectrum of 
lateral components of turbulence in the atmosphere over the frequency 
range affecting the dynamic response of an airplane: 


where k 1 = — and L is the integral scale of turbulence. The 
U 

value 2 jtL may be considered an approximate measure of the average 
eddy size in the turbulence. The data of reference 18 further indicate 
that L is on the order of 1,000 to 2,000 feet, and probably closer 
to 1,000 feet. 

The foregoing expressions, having asymptotic logarithmic slopes 
of -2 for values of k f <», were first suggested in reference 19 on 
the basis of wind-tunnel data and have since become generally accepted 
as reasonable expressions for defining the power spectra of lateral and 
longitudinal components of isotropic turbulence. Therefore, these 
expressions normalized by the mean-square values are used herein as the 
basis for calculations involving the components of atmospheric gusts as 
sensed at the center of gravity (point spectra) : 
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Its horizontal counterpart has the spectrum 
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In the calculations of reference 12, which involved the gusts sensed 
hy a wing having finite span, equations (46) and (47) were likewise used 
so that the data of reference 12 and this paper are in accord as to the 
spectra of gusts. However, equation (47) , as such, is not used directly 
in this paper and for this reason only equation (46) is plotted herein. 
The variation of equation (46) with k' is shown in figure 2. 


The equations which define the power spectral response of an air- 
plane to atmospheric turbulence (eqs. (18) to (21)) have been applied 
to three airplanes of varying size and wing span. The flight conditions 
and stability derivatives of the example airplanes are given in table I, 
and some physical dimensions necessary to the calculation of frequency- 
dependent force and moment coefficients are given in table II. Some of 
these physical dimensions require reference 11 for their concise defini- 
tion (sq, s-^, Xq, x^, and Xg) . 

By using the data presented in table I, the transfer functions of 
equations (26) to (35) were calculated. For discretely chosen values of 
reduced frequency a>', the in-phase and out-of-phase (real and imaginary 
parts of the transfer function) components were evaluated with a relay 
computer. Care was taken to choose as one point the natural frequency 
of the Dutch roll mode. 

By using the data of table II the complex moment and force coef- 
ficients (eqs. (38) to (45)) were evaluated at the same values of fre- 
quency as were used for evaluating the transfer functions. The rela- 
tionship 


was used to relate circular frequency in radians per second to reduced 
frequency in radians per span. The evaluation of certain of these coef- 
ficients requires further explanation. 


frequency are plotted in the figures of reference 11 for the example 
airplanes of this paper with one adjustment: 


APPLICATION OF METHOD 


Lateral Response of Three Airplanes 


The variations of 



C n per unit side gust (ref. 11) = U 
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The designation of example airplanes for this paper and reference 11 is 
consistent. 


The two rolling- and two yawing -moment coefficients due to Ug 
and w„ components acting on the wing were taken directly from the 

o 

figures of reference 12. Since these four complex moment coefficients 
represent the spanwise integration of the power spectra of the horizontal 
and vertical gust components, they appear in reference 12 only in the 
form of the absolute value squared, whereas equations (19) and (21) 
appear to require their real and imaginary components. However, by 
substituting equation (42) into equation (19 ) > the expression for the 
motion due to horizontal gusts becomes 



and since 


z*z = |z 


the equation becomes 
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where only the square of the absolute value of — appears and the 

\ u g/W 

product 
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g 


(due to horizontal gusts) 

W 8 1 


is plotted in figure 9 of reference 12. 


In the parallel evaluation of the motion due to vertical gusts, the 
substitution of equation (43) into equation (21) yields 
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where the product 


£l 

w„ 


K 

w g 



(due to vertical gusts) 


is plotted in figure 7 of reference 12 . 

Inasmuch as the figures of reference 12 are plotted against k', 
the relationship 


k . = ojL = cob = a/ ( 50 ) 

U p'U 3' 


where 



was used in the compatible choice of frequencies. Since the ratio of 
wing spans between the three example airplanes is on the order of 

the values of 3 ' were taken as 

3' A = 0.03125 

3' b = 0.0625 
3' c = 0.103 


which correspond to a value of L on the order of 1,100 feet. A 
rectangular span load distribution was used for airplane A and elliptical 
span load distributions were used for airplanes B and C. This choice was 
a matter of convenience since the plots of reference 12 indicate only a 
minor effect due to span loading distribution. 

The power spectral responses of the airplanes were completely cal- 
culated in three degrees of freedom, 0, t, and 3> "by using the power 
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spectra of the gust components given by equation (47) . Angular displace- 
ments of the airplanes are shown in figures 3 to 11 as follows: 


Airplane 

Angular displacement 

Figure number 



3 

A 


4 


3 

5 



6 

B 

* 
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3 

8 



9 

C 

♦ 

10 


3 

11 


In the (a) parts of figures 3 to 11 are plotted the three components 
and the sum for that particular angular displacement as expressed by 
equation (l8) . As may be seen from these figures, in none of the cases 
considered does the u component of gust velocity contribute perceptibly 
to the resulting motion. Therefore, in the (b) parts of figures 3 to H 
are plotted the relative contributions to the motion of the terms due 
only to the v and w components of the gusts. The (b) part of each 
figure then is a plot of equations (20) and (49) (where eq. (49) is a 
modified form of eq. (21)) and, as such, does not include any power- 
spectrum terms. For instance, the three components due to vertical 
gusts Wg shown in each (b) plot have been summed, multiplied by the 
power spectrum defined by equation (39) > and plotted in the (a) part 
of each figure. The u g and v g contributions were treated likewise 

with the exception that the components due to u g are not plotted in 
the (b) part of each figure. 

Figures 3 to 11 are plotted in this manner to permit observation 
of the end results without undue complication. When the relative magni- 
tudes of the contributing forces and moments are of particular interest, 
the (b) parts of the figures may be studied more closely. 
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Since in order to obtain the (a) parts of these figures, the 

|C I 2 

Wrr components are multiplied by — and the v K components 

s l w glW w g & 

multiplied by $ v , it should be noted that in the (b) parts of these 
& 

figures the units of the Wg components^are dimensionless while the 
v g components have the units * 

The 3 response plotted in figures 8, and 11 is the angular 
displacement with respect to the general air mass (that is, still air 
or fixed axes) and as such does not include the displacement with respect 
to the local air mass (gusts). 

From these example cases and figures , conclusions may be drawn on 
characteristics of method, trends, and further simplification of the 
calculations and a comparison of results from other methods of analysis 
may be made. 


Analysis 

Characteristics of the method .- In applying the method of this 
paper, several characteristics appear to be worthy of special attention. 
The first of these is the importance of calculating the response of the 
airplane at the exact frequency where the Dutch roll mode has a maximum 
amplitude. Since the Dutch roll mode is generally lightly damped and 
the square of the amplitude is required in the calculations, the power 
spectral response of the angular displacements of the airplane will 
generally exhibit a sharp peak in this frequency range. Enough datum 
points should be calculated in this region to insure an accurate repre- 
sentation of the response. 

A second characteristic of the method is the inaccuracy of equa- 
tion (20) as cd — » 0 for the $ and \|r modes of motion. The reason 
for this inaccuracy becomes apparent when the equation relating 0 
(or \|c) to v g 


v g 
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is written in the form 
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0 C Y/ 


The first half of equation (52) now represents the steady-state response 
to a uniform change in side gust over the airplane. The second half 
represents the additional response due to the nonuniform distribution 
of the gusts over the airplane. If the frequency-dependent expressions 

for the transfer functions i-, and as given by equations (2 6) , 

G n g Y 

(29), and (32), respectively, are substituted into equation (52) and the 


coefficients having the same powers of 

00 


0) 




are grouped ^£(a>) may be 
expanded in the form ^ a m aPj , the coefficient of the lowest order power 

m=l ' 

of co in the numerator may be shown to be identically equal to zero. 

In the numerical evaluations at low frequencies where the lowest order 
terms of equation (52) (or the square of the absolute value of eq. (52)) 
are predominant, small differences between large numbers appear unless % 
these terms are canceled beforehand. Since in the present method this 
cancellation is not convenient, a scattering of calculated values for 

and — begins to occur at some frequency below that of the Dutch 

V S V S CS Mr 

roll mode. However, since this scattering occurs only as — and — 

v g v g 

approach a constant (at go = 0) , this constant may he evaluated and the 
curve faired to that value. In figures ^y 6, 7> 9 y and 10 this 
fairing is denoted with dashed lines. In the range of frequencies where 
the response in (f and \| r due to Vg is faired, the total response of 

the three airplanes is due almost entirely to the Wg component of the 
gusts. 


Effect of trim angle of attack .- In analyzing the motions of the 
three example airplanes, the effect of small differences in trim -lift 
coefficients or trim angles of attack appears to he important. Although 
the trim angles of attack for all three airplanes are low, their rela- 
tive magnitudes are 


B:A:C = 1:1.36:2.50 
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Thus, airplane C was at an angle of attack from two to two and a half 
times as great as airplanes A and B. The stability derivatives for the 
wing alone in incompressible flow are functions of angle of attack (and, 
hence, lift coefficient) of the order (see ref. 20) : 

Cjp = Constant 

Cnp « cl 

C 7 cc a 
lr 

2 

Cn r * a 311(1 C Do 

The effect of increased angle of attack on the ratio of yawing moment to 
rolling moment may he seen to increase by 

c n p c n r Cdo 

— — oc a a a and 

C 7 C 7 00 

6p 

These ratios may explain trends such as that shown in the case of air- 
plane C for which the response in and (3 to w g at the higher 
frequencies is primarily due to the ability of Wg to produce yawing 

moment, whereas for airplanes A and B the response is primarily due to 
the ability of Wg to produce rolling moment. 

Effect of horizontal gusts .- While there is, in general, some cross 
power (or cross correlation) between any two velocity components measured 
at any two points in the turbulence, only those components lying in the 
XY-plane of the airplane are considered herein since the Z-gradients of 
the gusts are unimportant to the airplane. The airplane will, however, 
sense the correlation between the u g component acting at one point on 
the wing (at any instant) and the Vg component acting at some other 
point on the wing or vertical tail (a vortex effect) . It is the addi- 
tional effect of this correlation that has been neglected in the method 
of this paper. While this correlation is probably small, the fact that 
Ug has a negligible effect on the airplanes considered is reason enough 
for neglecting any cross -power effects involving Ug. In cases where 
the Ug component of gusts does have a pronounced effect on the airplane, 
the contribution of the cross power must be determined or further justi- 
fication be made for neglecting it. 

Response characteristics of example airplanes .- The calculated 
lateral responses of all three example airplanes exhibited a number of 
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characteristics in common. Most important of these was the total motion 
in each of the three degrees of freedom. 

In each case the response in roll was due primarily to Wg at low 
frequencies, Vg near the Dutch roll frequency, and Wg at the higher 
frequencies. The response of the example airplanes in yaw was due pri- 
marily to Wg at very low frequencies and due primarily to Vg at the 
Dutch roll and higher frequencies. The response in 3 of the three 
example airplanes was due primarily to v throughout the frequency 

o 

spectrum. 


Simplified Equations 


Based on the three conventional example airplanes, certain simplifi- 
cations to the response equations (l8) to (21) appear to he justified. 

For the low trim lift coefficients (or angles of attack) investigated, 
the effect of the horizontal gust (ug) was not discernible in any of 
the three modes of motion ((a) parts of figs. 3 to H) and could have 
been neglected. Furthermore, in every case investigated, the effect of 


g 


side force due to side gusts 
to 11) apparently because of the magnitude of 

1 

C 


— -) was negligible ((b) parts of figs. 3 


the A, JL, 
Cv (X 


Cy and the magnitude of 


and — transfer functions for the example airplanes; 

-/y ^y Cy 

however, these orders of magnitude are typical of most present-day 
airplanes. Other parameters appear to be negligible in some cases. When 
computing the response of the airplane in sideslip, neglecting the effect 
of Wg as well as that of u g appears to be justified. Although for some 

cases, for example, airplanes A and B, the effect of — — appears to be of 


w 


g 


small order in all modes, this effect is not generally true, as is shown 
by airplane C. As pointed out in a preceding section, the response in 
sideslip was due primarily to v g . 


In general, then, the preceding discussion leads to the simplifi- 
cation of equations (l8) to (21) as follows, where 0 may again be 
replaced by any of the lateral degrees of freedom in the form of angular 
displacements, velocities, or accelerations: 


= 


0 

2 

K + 

0 

V g 

v g 

V g 


■w 


g 


2 


(53) 


34 


NACA TN 3954 


and 
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COMPARISON OF METHODS 


Both airplanes A and B have Been analyzed in other papers (refs. 6 
to 9 ) at the same flight condition as that used herein for their responses 
to atmospheric turbulence by using methods which treat gusts as equiva- 
lent motions of the airplane in still air. These gusts are denoted as 

chr K 

side gusts 3 g , rolling gusts D0 g or -^- 2 , and yawing gusts Di|r or 

dug s 

-r— . The analysis of airplane A given in reference 6 is based on the 

oy 

power spectra derived from the analysis of reference 9 . Although the 
analysis of reference 6 also incorporates a refinement by considering 
the lag of the vertical tail penetrating the gust, the effect proved to be 
negligible on such a small fast airplane. Since in this analysis (ref. 6 ) 
the rolling gust affected only the wing, a comparison with the analysis 
of this paper (where the vertical gust likewise affects only the wing) 
may be made. Therefore, a comparison is made herein by using airplane A 
to illustrate any differences in the results of the two methods. 


Comparison of Gust Power Spectra 

Since the power spectrum of the v g component of atmospheric turbu- 
lence as used in references 9 and 6 and that used herein have in common 

O 

an asymptotic value of k" at large values (relatively speaking) of k, 
the arbitrary constant A, given in references 9 and 6 , is assigned a 
value such that the referenced v g spectrum and the v g spectrum of 
this report are asymptotically congruent. With this same constant, the 
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Sw 

D (f)g and spectra are compared, and all comparisons of the calcu- 

oy 

lated motions of airplane A are likewise based on this same constant. 


Side gusts .- The power spectrum of side gusts expressed by equa- 
tion (47) written in the form 


*V. 


g _ 


®P, 


g 


g 


Vy 


u 


_L_ 

«U 


1 + 


i~) 

\u 1 


i + == 


oiL^ 
U / 


(56) 


is comparable to that shown in figure U of reference 6 except for a 
constant. That is. 





(57) 


where 



(58) 


Plotted in figure 12 are the side-gust spectrum given by equation (56) 
and that of reference 6 with the constant K of equation (58) given the 
value 


K = 


72L 


(59) 


in order to get the agreement shown. Since both spectra are multiplied 
by the factor L, this agreement is independent of the choice of magni- 
tude for L. 

Rolling gusts .- Similarly, the power spectrum 
used in reference 6 may be compared with the power 
moment due to vertical gusts as given in reference 
Since 

Cl = - C Zp | = - ^ DSZS (60) 


of rolling gusts D^g 
spectrum of rolling 
12 and used herein. 


then 


3 6 
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and for airplane A, under the conditions considered herein, 

E(w g ) 


WcPlG - 2 


"g JJ ''Zp 
g Cz p 2 rt U : 


where E(wg) is the curve of 


2 tj3 
<J€ z (k') 


w, 


g 


^i p 2 /u 3 * 


for 3' = 0.03125 from 


reference 12, figure 7(a). If %(f, & is expressed with the same normali- 
zation used for the side gust (see eq. (56)), then 


g _ 4L 


w 


g 


= E(w g ) 

2 /y2 XU ^ 


(62) 


Equation (62) is plotted in figure 12 together with the power 
spectrum of D0 g given in reference 6 again with 


(63) 

rt u 2 

Since the spectra of reference 6 are the same as those given in refer- 
ence 9 except for a slight difference in frequency range, no difference 
is shown in the plots except that the range of data of reference 9 is 
denoted with ticks. 

Likewise shown in figure 12 is the spectrum for rolling gusts 
as predicted by Decaulne in reference 8. The theory is based on the 
assumption that the variation of vertical gusts across the span (dw g /ay) 
can be approximated with a constant gradient (a linear variation) and has 
as a spectrum r- -|2 

cob 


dWg 

6y 


which for small values of reduces to 



aib 

36 

sin 2u 

b 2 

(cn b\2 
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cos 
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^wg(-) 

dy 
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Based on this theory, when cu is less than the frequency corre- 
sponding to a wavelength of L (a value of L of 1,128 feet used in 
fig. 12), the vertical-gust spectrum will become a constant and the 

dwg/dy spectrum will fall off to zero at the rate of a) as a) 
approaches zero. This result is denoted by dashed lines in figure 12. 

At the higher frequencies shown in figure 12, the calculated spec- 
trum of reference 8 predicts almost no attenuation of rolling power of 
the wing with frequency. The spectrum of the same quantity according 
to references 9 and 6 would predict an attenuation for the wing of air- 
plane A, the attenuation beginning at about 8 radians /sec with the power 
falling off rapidly thereafter. The more comprehensive theory of refer- 
ence 12 shows a more pronounced attenuation of power over the plotted 
frequency range although all theories have approximately the same value 
at the frequency corresponding to the value of L. 

On the basis of this comparison, the concept of constant gust gra- 
dients across the span would appear to be insufficient, despite the 

agreement in the region of the break frequency The more com- 

prehensive theory of random gusts across the span predicts no attenuation 
below the break frequency such as is predicted by the theories of refer- 
ences 9 and 6. At frequencies greater than the break frequency, the 
attenuation predicted by the two theories again becomes inconsistent. 


Comparison of Lateral Response Power Spectra 

With the assigned value for the constant A given by equation (63) , 
the responses of airplane A in 0, \|r, and (3 0 due to 3 and D (f as 

o o 

calculated in reference 6 (fig. 10(a)) are plotted in figures 13(a), (b) , 
and (c) , respectively. Also replotted in figure 13 are the responses 
in (f), i|r, and 3 of airplane A due to v g and w g as obtained by the 

more comprehensive method of this paper (p 0 of ref. 6 is equivalent to 
3 of this paper) . Since reference 6 did not consider the Ug component 
of gust velocity, no comparison is made of the contribution of this 
component. Flight conditions are the same in both cases. 

In reproducing the frequency-response calculations of reference 6, 
the calculated transfer functions of \lr/D0g were found to be in error.; 

however, this error did not affect the conclusions made therein. The 
power spectral response of (^) D0 g shown in figure 13(b) has therefore 

been corrected and plotted for a somewhat larger frequency range than is 
given in reference 6. In the extended frequency range the D$ g spectrum 
which would be predicted by equation (65) is used. 


38 


NACA TN 395^ 


At the lower frequencies the differences in the response of the 
airplane due to side gusts as calculated by the two methods are due 
principally to the differences in the spectra of side gusts used by the 
two methods, as shown in figure 12. In the intermediate frequency range 
where the Dutch roll mode is predominant, a large difference in peak 
amplitudes is indicated by the two methods. This difference appears to 
be due to failure to evaluate the transfer functions at the exact fre- 
quency at which the Dutch roll mode has maximum frequency content 
(cpn = 3.16 radians /sec for ref. 6 as compared with ojq = 3-08 radians /sec 
for this paper) . In power spectral analyses where frequency amplitudes 
are squared, this factor can result in appreciable differences , as demon- 
strated in figure 13 . 

In the higher frequency range the differences in the theories are 
more pronounced since in this range the calculated effectiveness of the 
vertical gusts acting on the wing to produce lateral moments differs 
with the two theories. As pointed out in reference 11, the distribution 
of gust velocities along the fuselage begins to become an important factor 
in calculating the yawing moment and side force due to side gusts in this 
frequency range. 


CONCLUDING REMARKS 


By resolving the gust velocities of random isotropic atmospheric 
turbulence into three components alined with the stability axes of an 
airplane, the lateral forces and moments on the airplane and the resulting 
lateral motions of the airplane have been derived. The method involves 
the computation of the statistical forces and moments experienced by the 
airplane due to isotropic turbulence having random variations along the 
fuselage and across the wing span. On the basis of these gust-induced 
forces and moments referenced to the center of gravity of the airplane, 
the power spectral response of the airplane in any lateral degree of 
freedom (be it displacement, velocity, or acceleration) may be computed 
by using the appropriate transfer functions of the airplane. Since the 
gust velocities, the forces, and the moments are defined along the flight 
path in terms of their statistical average, the resulting motions are 
defined in terms of their power spectra, the mean-square value of the 
gust velocities, and the scale of atmospheric turbulence. 

By using this method, calculations have been made on three airplanes, 
which differed in size, experiencing continuous atmospheric turbulence 
while in a trimmed condition of low angle of attack. On the basis of 
these examples, certain possible simplifications to the equations have 
been brought out. Primary among these are the negligible effect of hori- 
zontal gusts Ug and the minor contribution of side force in computing 
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the motions. In general, then, the lateral motions at low trim angles 
of attack are due primarily to the ability of the side gusts Vg and 
vertical gusts w_ to produce yawing and rolling moments on the airplane. 

o 

With regard to lateral motions the present method of considering 
random variations of the gust velocities along the fuselage and across 
the span is a concept hitherto only approximated with constant gradients 
of gust velocity or completely ignored. Such approximations do not 
appear to he justified for gusts having wavelengths shorter than twice the 
span or fuselage length of the airplane. Although these less comprehen- 
sive methods correctly predict the trends in the vicinity of the Dutch 
roll mode of the airplane, they lead to the erroneous conclusion that 
the power spectra of the lateral motion due to the unsymmetrical compo- 
nents of the gust across the span (dwg/Sy and 8ug/8y) fall off to zero 
at wavelengths greater than the integral scale of the turbulence. The 
magnitude of difference between the methods increases with frequency 
above the Dutch roll mode since in this general region of frequencies 
the short-period (with respect to the size of the airplane) random varia- 
tions of gust velocity along the fuselage and across the span have a 
predominant effect on the moments and motions of the airplane. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January l8, 1957* 
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APPENDIX A 


PROOF OF POWER SPECTRAL AND CROSS-POWER 
SPECTRAL RELATIONSHIPS 


The physical and theoretical aspects of statistical methods of 
analyzing dynamic systems (generalized harmonic analysis) have been 
treated in a number of excellent papers. (Among these are refs. 13, 

14, and 21.) The following definitions and proofs are purely a mathe- 
matical means of utilizing the statistical concept and should be treated 
as useful supplements to the referenced statistical theory. 

Consider the simple equation 


9 = 



(Al) 


where 9, 6, and 7 are complex variables and 9/& and 9/7 are 

complex transfer functions such that 


9 - R 9 + i I 0 
5 = Rg + ilg 
7 = R7 + i I7 


9 

5 

9 

r 


,i0( 9/5) 


. 10 ( 8 / 7 ) 


(A2) 


With the use of an asterisk to denote the conjugate of a complex quan- 
tity, it is also true that 



Thus 
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By expansion in terms of equations (A2) it may be seen that 


&7 = 


/ * v* 

(6 7) 


0/0 

S\7. 


9\*e 

8/ 7 


and, in fact. 


9/9 


* * 


5 7 * ' 


in ^ 


(A5) 


Since the sum of a complex number and its conjugate is equal to twice the 
real part, that is , 

(R + il) + (R - il) = 2R 

then the substitution of equation (A5) into equation (Ah) yields 



2 

e 

2 


2 

0 

2 


2 

q \ *9 * 

0 

= 



6 

+ 



7 

+ 2R 

- - 6 7 



8 




7 




\8 / 7 


(A6) 


It is equally true to write 


2R 


0Y*0 * 

sir 57 


= 2R 


e/e\* * 
6 7 8r 


(A7) 


so that either form may he used. 

Dividing both sides of equation (A 6) by T and taking the limit 
of the statistical variables as T approaches infinity give the form 


iim s 

T— > 00 L 


lim — 

T^oo T 


lim ± 


T-» 00 


T 


+ 2R 


ifilim U*r 

5/ 7 T_*oo T 


(A8) 


where, by the definitions of equations ( 5 ) and (6) of the text, equa- 
tion (A8) is equivalent to 


$ 0 = 



2 


2 

" , v * 

9 

5 

% + 

CD| ^ 

$ r + 2R 

(t) 7 \ 


(A9) 
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If the two variables 7 and 6 are statistically independent, 
their cross power is zero. For equations with more than two variables 
the process is identical and, of course, yields more terms. For an 
equation whose solution is a function of m variables, there will 


result m power-spectrum terms and 


1 ^ 


cross-power terms in the 


n=l 

form of equation (A9) , or in the equivalent form 




2 


2 

/ \ * 

II 

CD 

e 

e 

5 

*8 + 

9 

7 

+ 2R 

•»] 


based on the identity of equation (A7) . 


(A10) 


For the cross spectrum between two statistical variables not 
appearing in the same equation, consider the equations 


0 



7 


and 

e = f 5 + ^ 7 (All) 

where the complex variable e, in addition to the definitions of equa- 
tion (A2) , is defined by H 


e 



+ il 


e 


Multiplying equation (A 3 ) by equation (All) gives 


(A12) 


9 € = ( -r 1 — 5 8 + - £ 


5/ 6 


7/ 7 


* / 0 \*e * / 9 \*€ * 

r7+ UJ7 6r+ 7 t 78 


(A13) 


where the identity 


(if; *•> ■ 


K;' v 


(A14) 


may be proved by the use of equations (A2) and (All) 


* 

9 € 




+ 2R 



, so that 



(A15) 
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Dividing both sides by T and taking the limit of the statistically 
dependent variables as T approaches infinity yield 


lira - 0 e = _ | lim ±|b 

T-»» T V 6 / 5 


'I 2 + (y)*r T|r 


+ 2R 


7) I h im A 

' ' T-*CO 


(A16) 


By the definition of equations (5) and (6) of the text, equation (Al6) 
becomes v 7 


eV*. 


9\*e 


•e« - (sj 8 % + (y) 7 ® r + 2R 

where again the identity 


f)*f V 


(A17) 






2R 

[(?) \ v] 

S 2R 

to n 


may be used in equation (A17) when desired. As before, $^5 = = 0 

where y and. 5 are statistically independent; but, in general, where 
he two equations are both functions of m variables, there will result 

m power -density terms and ^ (n-l) 

n=l 


cross -power terms. 
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APPENDIX B 

RELATIONSHIP BETWEEN POWER SPECTRA OF AIRPLANE FORCES AND 
MOMENTS AND POWER SPECTRA OF GUST VELOCITIES 


Consider the system defined by equation (2) of the text, which is 
rewritten here for convenience in the matrix form 



(Bl) 


where Cj is a force or moment coefficient (Cj, C n , or Cy) and aj m 
is a transfer function relating Cj to u^ The gust velocities Uy, 
u 2 , and u^ where m = 1, 2, or 3) are the three orthogonal velocity 

components of isotropic turbulence, and the cross power between any two 
components is zero; that is. 


u l u 2 


= <D. 


u 2 u 


= «, 

3 u 3 




= 0 


(B2) 


In the power- spectral domain (see appendix A) the first equation 
of the system is defined by 


*Ci = 


an 


'*11! + 


ai2 


*uo + 


ai3 


*u^ + 


2R 


a ll a 12 a> u 1 u 2 + a 12 a 13 $ u 2 u 5 + a 13 a ll 4 u 3 U! 


(B3) 


All the cross-power terms of equation (B3) are zero. The other equa- 
tions yield similar results; therefore, the system is defined by 



(B4) 


which is the result given by equation (9) of the text. 
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Now consider the cross-power terms of the Cj matrix, where, for 
instance. 


* 0^2 = 


lim J C;l*C 2 
T— >00 


= lim i 

T-*x> 


3 
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3 

'r~< 

/ a lm u m 

L — j 
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lim if \ a i in *a 2 rn n m *u m + 6 cross -power terms 


involving u-^ U 2 , U 2 u^, and so forth 


In the limit, then 


<I>C l C 2 = X aijn * a2m J im T + 0 


m=l 
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“ ^ a lm a 2m ®u m 


m=l 


where by definition 
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A similar derivation for the other two cross-power combinations yields 
similar results; therefore, in general. 


|®CjC k 


* 

a jm a km 


Um 


(jA) ( B8a ) 


where C^. is also a force or moment coefficient and a ^ is a transfer 
function relating to u m , and 


| 4 0jC k 


1 J m l 




(j=k) (B8b) 


Equation (B8b) is seen to be identical to equation (B4) . In expanded 
form, for the case considered 


1 C 2 


* 

a ll a 21 

* 

a~| p ^22 

* 

a 13 a 23 

! 2 C 3 
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a 21* a 31 

a 22 a 32 

* 
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* 
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* 
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r w 


g 


(B9) 


In the terms of equation (2) of the text, this result is equivalent to 
equation (10) of the text. Furthermore, other equally correct arrange- 
ments are possible, since 


(BlOa) 

(jA) (biod) 



, \ 

^ ®C k Cj| =' 

j*C k ' 


®CvC 


k u j 


P CkCj 


ajm^-km 




*u, 


L m 


Ufll 


(j=k) (BlOc) 


Equation (BlOc) is identical to equations (B4) and (B8b) . 
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TABLE I.- FLIGHT CONDITIONS AND STABILITY DERIVATIVES OF EXAMPLE AIRPLANES 


(a) Flight conditions 


Flight condition 

Airplane 

A 

B 

c 

v ^ 

50,000 

4,000 

55,000 

U, ft /sec 

69 6 

318 

700 

W, lb 

12,600 

28,000 

125,000 


50 

10.9 

51.85 


0.242 

0.177 

0.445 

tan 7 

0 

0 

0 

a Q , radians 

0.0586 

0.0525 

0.0825 


(b) Stability derivatives 


Stability derivative . 

Airplane 

A 

B 

C 

Total airplane 

K ^ 

K x 

0.01485 

O.OI58 

0.0511 

Kz 2 

0.0504 

0.0500 

0.072 

% 

-0.00062 

-O.OOI58 

0 

% • • 

-0.45 

-0.4875 

-0.44 

c i r 

0.040 

0.1055 

0.149 

% 

-0.11 

-0.1084 

-0.14 

s 

-0.055 

-0.0508 

-0.0275 

Cn r ’ 

-0.15 

-0.1041 

-0.156 

% 

0.12 

O.O851 

0.12 

Oy 

*P 

-0.015 

0 

0 

C Y r 

0.225 

0 

0 

% 

-0.58 

-0.765 

-0.6l 

c Ax 

4.15 

5.44 

5.38 

Wing 

n P / l V 

0.025 

0.0208 

0.054 

C n rl°lr 

-0.075 

-0.055 

-0.052 



-0.0475 

-0 . 0 666 

-0.0876 

Cdo 

0.01 

0.006 

0.01 

Vertical tail 

CQ 

-0.42 

-0.337 

-0.584 

do/&$ 

0.19 

0.242 

0.215 

% 

0.185 

0.129 

0.154 

C, 

-0.0625 

-0.0418 

-0.0524 
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TABLE II.- PHYSICAL DIMENSIONS 
OF EXAMPLE AIRPLANES 


Dimension 

Airplane 

A 

B 

c 

b, ft 

35-25 

70 

116 

S, sq. ft 

250 

540 

1,428 

Sip, sq. ft 

55 

71.35 

230 

A 

4.98 

9.07 

9.4 

T, deg 

4.0 

4.50 

3.6 

A . . 

0.46 

0.45 

0.42 

h z , ft 

4.4 

7.0 

13.0 

l t > ft 

14.8 

27.6 

46.6 

* 0 ; ft 

18 

15.25 

48.5 

x l> ft 

5-7 

22.0 

36.6 

x 2 , ft 

20 

29.4 

51.6 

s 0 , ft 

2.7 

3.0 

5-3 

s l> ft 

8.5 

10.6 

18.5 




u q 


(a) Turbulence structure and mean orientation 
of stability axes of airplane. 


VJ i 


Figure 1.- Sign convention and stability axes of airplane flying through 

atmospheric turbulence. 


+l£6£ Mi VOVM 


^earth 



(b) Inst ant aneous orientation of stability axes of airplane with respect 
to earth and body axes. Positive sense of axes , velocities, and angles 
is indicated by arrows. 


Figure 1.- Concluded. 
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Figure 2.- Power spectral density of u, v, and w components of homogeneous isotropic turbulence 

as sensed at single point in turbulence. 
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Normalized power spectral density of roll angle, <£, due to gusts, 
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Frequency, cu, radians/sec 

(a) Components and sum of components of power spectral density of roll 
angle due to three components of gust velocity . 

Figure 3.- Response in roll angle of airplane A flying through continuous 

atmospheric turbulence . 
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(b) Roll angle due to forces and moments induced by v g and w g components 
of gust velocity. Shown are sum and relative magnitudes of parts. 


Figure 5*- Concluded. 


Normalized power spectral density of yaw angle, i//, due to gusts, 
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Frequency, <0, radians/sec 

(a) Components and sum of components of power spectral density of yaw 
angle due to three components of gust velocity. 


Figure 4.- Response in yaw angle of airplane A flying through continuous 

atmospheric turbulence. 
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(b) Yaw angle due to forces and moments induced by v g and w g components 
of gust velocity. Shown are sum and relative magnitudes of parts. 

Figure 4.- Concluded. 
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(a) Components and sum of components of power spectral density of side- 
slip angle due to three components of gust velocity. 

Figure 5 .- Response in sideslip angle of airplane A flying through con- 
tinuous atmospheric turbulence. 
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(b) Sideslip angle due to forces and moments induced by Vg and 
Wg components of gust velocity. Shown are sum and relative 
magnitudes of parts. 


Figure 5«- Concluded. 
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(a) Components and sum of components of power spectral density of roll 
angle due to three components of gust velocity. 


Figure 6.- Response in roll angle of airplane B flying through contin- 
uous atmospheric turbulence. 
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(b) Roll angle due to forces and moments induced by v g and w g components 
of gust velocity . Shown are sum and relative magnitudes of parts . 


Figure 6.- Concluded. 


Normalized power spectral density of yaw angle, i//, due to gusts, 
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(a) Components and sum of components of power spectral density of yaw 
angle due to three components of gust velocity. 


Figure 7 - - Response in yaw angle of airplane B flying through continuous 

atmospheric turbulence . 
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(b) Yaw angle due to forces and moments induced by Vg and Wg components 
of gust velocity. Shown are sum and relative magnitudes of parts. 


Figure 7 *- Concluded. 
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(a) Components and sum of components of power spectral density of side- 
slip angle due to three components of gust velocity. 


Figure 8.- Response in sideslip angle of airplane B flying through 
continuous atmospheric turbulence. 
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(b) Sideslip angle due to forces and moments induced by v g and 
w components of gust velocity. Shown are sum and relative 

s 

magnitudes of parts . 

Figure 8.- Concluded. 
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(a) Components and sum of components of power spectral density of roll 
angle due to three components of gust velocity. 

Figure 9.- Response in roll angle of airplane C flying through 
continuous atmospheric turbulence . 
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(b) Roll angle due to forces and moments induced by Vg and 
Wg components of gust velocity. Shown are sum and relative 
magnitudes of parts. 


Figure 9.- Concluded. 


(radians) 2 /(ft/sec) 2 ft 

Normalized power spectral density of yaw angle, vj/, due to gusts, radians/sec 
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(a) Components and sum of components of power spectral density of yaw 
angle due to three components of gust velocity. 

Figure 10.- Response in yaw angle of airplane C flying through 
continuous atmospheric turbulence . 
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(b) Yaw angle due to forces and moments induced by Vg and 
Wg components of gust velocity. Shown are Siam and relative 
magnitudes of parts . 


Figure 10.- Concluded. 
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(a) Components and sum of components of power spectral density of side 
slip angle due to three components of gust velocity. 


Figure 11.- Response in sideslip angle of airplane C flying 
through continuous atmospheric turbulence . 
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(b) Sideslip angle due to forces and moments induced by Vg and 
Wg components of gust velocity. Shown are sum and relative 
magnitudes of parts. 


Figure 11.- Concluded. 


v 2 /L )2 w 2 /u 2 ' radians/sec 



Frequency, oj, radians/sec 


Figure 12.- Comparison of several theories for power spectra of side gusts and rolling gusts as 
experienced by airplane A. U = 696 ft/sec; b = 35*25 feet; L = 1,128 feet. 
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Frequency, <0, radian/sec 


(a) Power spectral response in roll angle </> of airplane A. 

Figure 13*- Comparison between the lateral responses of airplane A in 
atmospheric turbulence as determined by method of present paper and 
method of reference 6. 
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(b) Power spectral response in yaw angle \|r of airplane A. 
Figure 13.- Continued. 


radians/sec 


1Y 


NACA TN 395U 


75 



Frequency, ai, radians/sec 


(c) Power spectral response in sideslip angle p of airplane A. 

Figure 13 •- Concluded. 
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